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We present the first analytical inspiral-merger-ringdown gravitational waveforms from binary black holes 
(BBHs) with non-precessing spins that is based on a description of the late-inspiral, merger and ringdown in 
full general relativity (GR). By matching a post-Newtonian description of the inspiral to a set of numerical 
calculations performed in full GR, we obtain a waveform family with a conveniently small number of physical 
parameters. The physical content of these waveforms includes the "orbital hang-up" effect, and "spin flips". 
These waveforms will allow us to detect a larger parameter space of BBH coalescence (including a considerable 
fraction of precessing binaries in the comparable-mass regime), to explore various scientific questions related to 
gravitational- wave astronomy, and could significantly improve the expected detection rates. 



Introduction. — Coalescing black-hole (BH) binaries are 
among the most promising candidate sources for the first di- 
rect detection of gravitational waves (GWs). Such observa- 
tions will lead to precision tests of the strong-field predictions 
of general relativity as well as provide a wealth of information 
relevant to fundamental physics, astrophysics, and cosmology. 
Computation of the expected waveforms from such sources is 
a key goal in current research in gravitation. 

While the inspiral and ring-down stages of the BH coales- 
cence are well-modeled by perturbative techniques, an accu- 
rate description of the merger requires numerical solutions of 
Einstein's equations. Although performing numerical simu- 
lations densely sampling the entire parameter space of BH 
coalescence is computationally prohibitive, waveform tem- 
plates modeling all the three stages can now be constructed 
by combining analytical- and numerical-relativity results, dra- 
matically improving the sensitivity of searches for GWs from 
BH binaries, and the accuracy of estimating the source param- 
eters 03-2). To date, inspiral-merger-ringdown (IMR) tem- 
plates have been computed only for nonspinning BH bina- 
ries HI [3}|6l, which are employed in GW searches, and in a 
number of astrophysical studies. However, nonspinning BHs 
are expected to be astrophysically rare, and most BHs in na- 
ture may be highly spinning QO. This necessitates the in- 
clusion of spinning-binary waveforms in GW searches. But, 
spin adds six parameters (three components for each BH), and 
each additional parameter in a search template bank leads to 
a higher signal-to-noise-ratio (SNR) threshold for a confident 
detection (9). Also, this requires sufficiently accurate numer- 
ical simulations across this large parameter space, which are 
not yet available. Moreover, implementing a search covering 
the full spin parameter space has proven to be difficult. 

In this letter, we present an IMR waveform family modeling 
the dominant harmonic of binaries with non-precessing spins, 
i.e., spins (anti-)aligned with the orbital angular momentum. 
Aligned-spin binaries are an astrophysically interesting popu- 
lation as such systems are expected from isolated binary evo- 
lution and in gas-rich galactic mergers ifTUl [Til . Such sys- 



tems also exhibit interesting strong-gravity effects like the 
"orbital hang-up" iTTSl and "spin flips" [13]. We make use 
of the degeneracies in the physical parameters to parametrize 
our waveform family by only the total mass M = m\ + mi of 
the binary, the symmetric mass ratio rj = mim^jM 1 , and a 
single spin parameter % = (1 + 5)%\/2+ (1 — 5)^/2, where 
8 = (mi — mi)jM and Xi = Si/mf, Sj being the spin angu- 
lar momentum of the ith BH. The last feature is motivated 
by the observation (see e.g., Tl4\ ) that the leading spin-orbit- 
coupling term in post-Newtonian (PN) waveforms is domi- 
nated by this parameter. We also show that this waveform 
family is able to capture a significant fraction of precess- 
ing binaries in the comparable-mass regime, the reason being 
that the phase evolution is primarily governed by the spin- 
orbit coupling, which is determined by the component of the 
spins along the angular momentum. This provides an efficient 
and feasible way of searching for precessing binaries in the 
comparable-mass regime. 

Numerical simulations. — Binary BH (BBH) waveforms 
covering at least eight wave cycles before merger were pro- 
duced by solving the full Einstein equations numerically, as 
written in the "moving-puncture" 3+1 formulation |[T5l [T6l . 
The numerical solutions were calculated with the BAM ifTTl 
[TH, CCATIE 020 and LLAMA J20) codes. Initial momenta 
were chosen to give low-eccentricity inspiral, using either an 
extension of the method described in ED . or the quasicircular 
formula used in [22 1. GWs were extracted at R ex = 9QM with 
BAM, R ex = 160M with CCATIE and at future null infinity 
with LLAMA, using the procedures discussed in IfTTl [T9l l23l . 
In all simulations the GW amplitude is accurate to within at 
least 10% and the phase to within at least 1 radian over the 
duration of the simulation. Most of the waveforms employed 
in the construction of the analytical templates are significantly 
longer (12-22 cycles) and more accurate [24|. 

Seven sets of simulations were used in this paper: (1) 
Equal-mass binaries with equal, non-precessing spins Xi — 
±{0.25,0.5,0.75,0.85}, described in Il24ll25l (2) Unequal- 
mass, non-precessing spin binaries with q = m\/m2 = 
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FIG. 1: Phenomenological parameters \ff2 and computed from 
equal-spin hybrid waveforms (dots), and analytical fits given by 
Eq. |2j (surfaces). Test-mass limit is indicated by black traces, 77 
is the symmetric mass ratio, % the spin parameter, and M the total 
mass. 



{2,2.5,3} and %i = {±0.5,0.75}. (3) Nonspinning bi- 
naries with q = {1,1.5,2,2.5,3,3.5,4}. (4) Unequal- 
mass, unequal-spin binaries with q = {2,3} and {X\,X2) = 
(—0.75,0.75). (5) Equal-mass, unequal-spin binaries 
with Xi = ±{0.2,0.3,0.4,0.6}. (6) Equal-mass, precess- 
ing binaries with spin vectors (0.42,0,0.42), (0,0,0) and 
(0.15,0,0), (0,0,0). (7) Precessing binary with q = 3 and 
spins (0.75,0,0), (0,0,0). The simulation sets (l)-(4) and (7) 
were performed with the BAM code, set (5) with the C CAT IE 
code, and set (6) with the LLAMA code. The analytical wave- 
form family is constructed only employing the equal-spin sim- 
ulation sets (l)-(3), while sets (4)-(7) were used to test the 
efficacy of the template family to model the expected signals 
from more general spin/mass configurations. Two additional 
waveforms were used to test the template family: the equal- 
mass, non-spinning simulation performed by the Caltech- 
Cornell group (26), and the q = 1.25 precessing binary sim- 
ulation with spins (-0.09, 0.48, 0.35) , (-0.2, -0. 14, 0.32) by 
the RIT group 11271 . 

Constructing hybrid waveforms. — Following ITJ |5J, we 
produce a set of "hybrid waveforms" by matching PN and 
numerical-relativity (NR) waveforms in an overlapping time 
interval [t\ , £2]. These hybrids are assumed to be the target sig- 
nals that we want to detect. For the PN inspiral waveforms 
we choose the "TaylorTl" waveforms at 3.5PN ll28ll phase ac- 
curacy, considering spin terms up to 2.5PN |29l [30l . This 
is motivated by PN-NR comparisons of equal-mass spinning 
binaries, in which the accuracy of the TaylorTl approximant 
was found to be the most robust Il24ll25ll . We include the 3PN 
amplitude corrections to the dominant quadrupole mode OTI 
and the 2PN spin-dependent corrections lfT4ll . which greatly 
improved the agreement between PN and NR waveforms. 
For precessing waveforms, the spin and angular momenta are 
evolved according to Il30ll32l . 

If h(f) = /i+(f) — ih x (t) denotes the time-domain wave- 
form from a binary, we match the PN and NR waveforms 
by doing a least-square fit over time- and phase shifts be- 
tween the waveforms, and a scale factor a that reduces the 
PN-NR amplitude difference. The NR waveforms are com- 
bined with the "best-matched" PN waveforms in the follow- 
ing way: h h 5*(f) = ai(t) h NR (f) + (1 - %{t)) h PN (f), where T 
ranges linearly from zero to one for t € [t\ , to]. 



Waveform templates for non-precessing binaries. — The 
analytical waveforms that we construct can be written in the 
Fourier domain as h(f) = A{f) e~ l ^^\ where 

C f-y 6 (l+zl 2 a l v') if/</! 

A(f) EE Vf- V6 \ W m /'- 2 / 3 (l+L2 =l£i y) if/ 1 </</ 2 

{ w r J?(f,f 2 ,o) if/ 2 </</ 3 , 
¥(/) ee 2nft + w + —^(l + Y,v k Vk ). (1) 

l/OTJV k=2 

In the above expressions, /' ee f/fi, v = (ttM/) 1 / 3 , £1 = 
1.4547£- 1.8897, e 2 = - 1.8 153 % + 1.6557 (estimated from 
hybrid waveforms), ^ is a numerical constant whose value 
depends on the sky-location, orientation and the masses, (X2 = 
-323/224 + 45177/168 and a 3 = (27/8 - 11 r}/6)% are the 
PN corrections to the Fourier domain amplitude of the (I — 
2,m = ±2 mode) PN waveform lfT4l . to is the time of ar- 
rival of the signal at the detector and <po the corresponding 
phase, j£?(/,/2, c) a Lorentzian function with width a cen- 
tered around the frequency fi, w m and w r are normalization 
constants chosen so as to make A(f) continuous across the 
"transition" frequencies fi and f\, and fa is a convenient cut- 
off frequency such that the power of the signal above this fre- 
quency is negligible. The phenomenological parameters % 
and fik ee {/1 ,/2, cj,/3} are written in terms of the physical 
parameters of the binary as: 

3 N 3 N v ('7) n i v j 1 „0 

i=l ;=0 i=l 7=0 llM 

where = min(3 — i,2) while x[ !; ' and are tabulated 
in Table [j] Figure [TJ plots an example of this map from the 
phenomenological- to physical- parameter space. 

We match these waveforms to 2PN accurate adiabatic in- 
spiral waveforms in the test-mass limit. In the 77 — > limit, 
the phenomenological parameters reduce to the following: 

/W/lso- V*->¥#. ( 3 ) 

where f^ so and /q NM are the frequencies of the last stable or- 
bit ll33ll and the dominant quasi-normal mode, and Q° is the 
ring-down quality factor |34| of a Kerr BH with mass M and 
spin x, while ly® are the (2PN) Fourier domain phasing coef- 
ficients of a test-particle inspiralling into the Kerr BH, com- 
puted using the stationary -phase approximation 1 14 1. 

The test-particle-limit waveforms suffer from two limita- 
tions: 1) we assume that the evolution of the GW phase at 
the merger and ringdown stages is a continuation of the adia- 
batic inspiral phase, and 2) in the absence of a reliable plunge 
model, we approximate the amplitude of the plunge with 
y/-2/3 ^ +Y^=i £;v ! ). Nevertheless, in the test-mass limit, it 
is expected that the signal will be dominated by the long inspi- 
ral stage (followed by a quick plunge and ringdown), and the 
inspiral is guaranteed to be well-modelled by our waveforms. 
More importantly, the imposition of the appropriate test-mass 
limit in our fitting procedure ensures that the waveforms are 
well behaved even outside the parameter range where current 
NR data are available. Because of this, and the inclusion of 
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FIG. 2: Top and right plots: Match and FF of our analytical IMR 
templates with equal-spin hybrid waveforms constructed from simu- 
lation sets (l)-(3). Bottom left: FF of non-spinning IMR templates 
proposed in 1 1 6] with the equal-spin hybrid waveforms. A compar- 
ison with the other plots demonstrates the effect of neglecting spins. 
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FIG. 3: Top panel: Match and FF of the our templates with unequal- 
spin hybrid waveforms constructed from simulation sets (4) and (5), 
and the Caltech-Cornell non-spinning simulation. Bottom left: FF 
with precessing hybrids constructed from sets (6) and (7), and the 
RIT simulation. Bottom right: Fraction of generic precessing PN 
waveforms (M = 20M a ) producing fitting factor FF with the tem- 
plates proposed in this paper — 100% (85%) 65% of the binaries 
with q = 1(4)10 produce FF > 0.965. 



the PN amplitude corrections, these waveforms are expected 
to be closer to the actual signals than the templates proposed 
in ID |6l in the non-spinning limit. However, since the pa- 
rameter space covered by the NR simulations is limited, we 
recommend that these waveforms be used only in the regime 
q < 10 and —0.85 < % ^ 0.85. Also, these are meant to model 
only the late-inspiral, merger and ring down (M/gw > 10~ 3 ), 
i.e., signals from binaries in the mass-range where merger- 
ringdown also contribute to the SNR, apart from inspiral. 

We have examined the "faithfulness" P5l of the new tem- 
plates in reproducing the hybrid waveforms by computing the 
match (noise-weighted inner product) with the hybrids. Loss 
of the SNR due to the "mismatch" between the template and 



the true signal is determined by the match maximized over 
the whole template bank - called fitting factor (FF). The stan- 
dard criteria for templates used in searches is that FF > 0.965, 
which corresponds to a loss of no more than 10% of signals. 

Match and FF of the analytical waveforms with the equal- 
unequal-) spin hybrid waveforms are plotted in Fig.[2](Fig.[3]l, 
using the Initial LIGO design noise spectrum [36|. Note that 
the analytical waveform family is constructed employing only 
the equal-spin hybrid waveforms (Fig.|2j. The PN-NR match- 
ing region used to construct the unequal-spin hybrids (Fig. [3} 
are also different from that used for equal-spin hybrids. These 
figures demonstrate the efficacy of the analytical templates 
in reproducing the target waveforms - templates are "faith- 
ful" (match > 0.965) either when the masses or the spins 
are equal, while they are always "effectual" [35] in detection 
(FF > 0.965) E). In contrast, the bottom left plot of Fig. [2] 
shows the FF of the non-spinning IMR template family pro- 
posed in HI [6) with the equal-spin hybrid waveforms. FFs 
as low as 0.8 suggest that up to 50% binaries may go unde- 
tected if nonspinning IMR templates are employed to search 
for binaries with high (aligned) spins. 

Bottom left plot of Fig. [3] shows the FF and match of the 
template family with four precessing hybrid waveforms. The 
high FFs are indicative of the effectualness of the templates in 
detecting precessing binaries. Since presently not enough NR 
simulations are available to make a quantitative statement, and 
since we expect the effect of precession will be predominant 
in the case of lower mass binaries (when large number of cy- 
cles are present in the detectors band), we might be able to ac- 
quire some useful indication by studying precessing PN wave- 
forms. We have performed a Monte-Carlo simulation where 
we generate generic precessing "restricted" 3.5PN waveforms 
(with 2.5PN spin corrections) with M = 2OM , q = {1,4, 10} 
and spins uniformly distributed in a sphere of radius 0.98, and 
compute the FF with the templates proposed in this paper. 
Inclination of the binary's plane with the line of sight from 
the observer is also randomly selected from [0, |]. Bottom 
right plot of Fig. [3] shows the cumulative distribution of the 
FF, strongly indicating the effectualness of the templates in 
detecting precessing binaries in the comparable-mass regime. 
Almost all binaries producing FF < 0.965 have (the astrophys- 
ically less likely [ 10]) tilt angles > | for the massive BH, and 
undergo flipping of the orbital plane. 

Distance to optimally oriented BBHs producing optimal 
SNR of 8 at Initial LIGO noise spectrum is shown in Fig. [4] 
which demonstrates the dramatic effect of spin for detection 
of high-mass binaries; if most BBHs are highly spinning, then 
LIGO will be able to detect BH coalescences up to lGpc, thus 
increasing the event rates as much as five times compared to 
predictions based on models of nonspinning binaries. For Ad- 
vanced LIGO, the distance reach is as high as 20 Gpc. 

Conclusions. — We combine state-of-the-art results from 
analytical- and numerical relativity to construct for the first 
time a family of analytical IMR waveforms for BBHs with 
non-precessing spins from "first principles". These wave- 
forms are also able to detect a significant fraction of the pre- 
cessing binaries in the comparable-mass regime, and we show 
that for the purposes of GW detection it is sufficient to rep- 
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TABLE I: Phenomenological parameters describing the analytical waveforms (see Eq. l|2j). In test-mass limit, they reduce to the appropriate 
quantities given by perturbative calculations 1 14 33 34 1 . The test-mass limit of f\ is a fit to the frequency at the last stable orbit given in 1 33 1. 
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FIG. 4: Distance to optimally-located and oriented- equal-mass bi- 
naries with (equal) spin x producing optimal SNR 8 in Initial LIGO. 

resent the spins by a single parameter. This will considerably 
simplify the use of our waveforms in GW searches in the near 
future, and will accelerate the incorporation of NR results into 
the current effort for the first detection of GWs. There are 
many other immediate applications of our waveforms: injec- 
tions into detector data will help to put more realistic upper 
limits on the rate of BBH coalescences, and to compare the 
different algorithms employed in the search for BBHs, while 
employing these in population-synthesis studies will provide 
more accurate coalescence rates observable by the current and 
future detectors. Our method can readily be generalized to in- 
corporate non-quadrupole harmonics, larger portions of the 
BBH parameter space and further information from analytical 
approximation methods or numerical simulations. 
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